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Abstract 

A simple argument is presented which, based on the minimal coupling Lagrangian 
for a many-body system, keeps the gauge invariance of the many-body Schrodinger 
equation with explicit Coulomb potential. The elimination of longitudinal electric 
field does not necessarily lead to the breakdown of gauge invariance. The total time 
derivative term in the matter-EM field interaction in the Lagrangian is canceled 
out by the choice of Coulomb gauge. The remaining interaction is described by 
transverse vector potential and the longitudinal electric held which is the homo¬ 
geneous solution of Gauss law. This leads directly to the gauge invariant forms 
of the linear and nonlinear constitutive equations. It is discussed how to recon¬ 
cile this result with the cases of an isolated matter interacting with external charges. 


1 Introduction 

Quantum mechanical description of electromagnetic (EM) response of matter re¬ 
quires the use of scalar and vector potentials for EM held. A general formulation 
starts usually from the so-called minimal coupling Lagrangian C for the system of 
charged particles and EM held. It is the sum of particles Lagrangian and that of 
free EM held. The particle part is a sum of single particle Lagrangian in a given 
EM held, n^/2m — ecf) + {elc)v ■ A, over all the particles in consideration. (See 
eq.([T]) for notations.) The minimum action principle leads to [a] Newton equation 
of motion of each particle under Lorentz force, and [b] Maxwell equations for EM 
held. This fact guarantees the soundness of the Lagrangian. A gauge transforma¬ 
tion of the Lagrangian leads to an addition of a total time derivative to the original 
Lagrangian. Since the minimum action principle is not ahected by such a term, we 
get a gauge invariant EM response. 

The Hamiltonian corresponding to C is derived via standard procedure, and the 
Schrodinger equation in terms of this Hamiltonian allows us to calculate the expec¬ 
tation values of various physical quantities. Among them, that of induced current 
density gives the constitutive equation, which, together with Maxwell equations, 
makes up the fundamental equations to calculate the EM response of matter. In 
semiclassical theory where EM held is treated as classical (non-quantized) quan¬ 
tity, these fundamental equations are solved as simultaneous equations for a given 
initial condition of matter and EM held. Maxwell equations give EM held as a 
functional of current and charge densties, while constitutive equation gives current 
density as a functional of EM held. 

According to our standard conception about matter, the Coulomb potential 
among charged particles is a part of matter Hamiltonian. The energy levels of a 



matter are given as the eigenvalues of the matter Hamiltonian consisting of the 
kinetic energies and the Coulomb potential of all the particles in the system, as 
typically seen in the Rydberg series of a hydrogen atom. However, £ does not 
contain the Coulomb potential explicitly. It emerges from rewriting the self-energy 
of longitudinal (L) electric held in terms of Gauss law V ■ E = p/e^. When we 
derive the matter Hamiltonian with explicit Coulomb potential from the minimal 
coupling Lagrangian, the gauge invariance of the many-body Schrodinger equation, 
as mentioned in eqs. fl26|) . (HT} below, is not automatically kept. There seems to 
be no careful argument of this problem in literatures, to the author’s knowledge. 
Some people argue that gauge invariance is broken by eliminating the L electric 
held to derive Coulomb potential. In this note, a simple argument is given to keep 
the gauge invariance in the many-body Schrodinger equation. 

An additional point requiring a careful consideration is how we extract, from C 
for a global sysytem, the Lagrangian to describe the quantum mechanical motion of 
the isolated matter system (charged particles) of our interest. For that purpose, we 
need to write C in terms of the internal and external variables, and after rewriting it 
into Hamiltonian with explicit Coulomb potential, we drop the irrelevant variables 
of the external origin. Thereby, it is important to note that the L electric held in 
C contains, in adition to the chrage induced component, the solution of the 
homogeneous equation V ■ E = 0. This term plays an essential role for the gauge 
invariance of the many-body Schrodinger equation. 

Once we derive the matter Schrodinger equation in a gauge invariant form, the 
expectation value of an arbitrary physical quantity should be gauge-independent. 
Among all. Coulomb gauge is a special one in the sense that the Hamiltonian is 
already written in terms of the gauge independent components of EM held, i.e., 
the transverse (T) vector potential and homogeneous L electric held Eq^\ 
Hence the result obtained in Coulomb gauge directly gives the gauge independent 
form. This conclusion applies, not only to linear, but also to all the nonlinear EM 
responses. 


2 Formulation 

We start from the minimal coupling Lagarngian to describe a system of interacting 
charged particles and EM held as 

£ = ^ - ei(t){ri) + avi ■ A{re) } 

i ^ 

+|/dr , (1) 

where E = —dA/dt — V0, B = V x A, and mi,e£,ri,ve are the mass, charge, 
coordinate, and velocity, respectively, of the Gth particle, and 0 and A scalar and 
vector potential, respectively. The interaction terms in £ can be rewritten as 

+ CiVi ■ A{ri) }= f dr {-p(r)0(r) + J{r) ■ A(r)} , (2) 

e 

where 

p(r) = Y,ee5{r-re), (3) 
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are charge and current densities, respectively. 

As well known, the Lagrange equations due to this C lead to the Newton equa¬ 
tion of motion of the Ath particle 
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and the Maxwell equations for EM held 
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Maxwell equations give EM held as a functional of p and J, which also contains 
the solution of the homogeneous equations (for p = 0, J = 0). The existence of the 
held component independent of charge density is important for the gauge invariant 
rewriting of C into many-body Hamitonian with explicit Coulomb potential. 

Now we rewrite C by using Gauss law, V ■ E^^"^ = p/eo, which leads to Coulomb 
potential and eliminates the time derivative of from C. The general solution 
of Gauss law is the sum of the charge induced held S and the solution E^^^ of the 
homogeneous equation as 
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where Eq^^ is a spatially constant vector satisfying V ■ E^^^ =0. Its representation 
in terms of scalar and vector potentials E^^^ = —(dA^^^/dt) — V0o has inhnite 
choices according to gauge transformation. 

Since we are interested in an isolated matter system, we divide p into the internal 
and external components 

P = Pint + Pext (10) 

where pi^t belongs to the isolated system. Correspondingly, the charge induced 
held is also the sum of the components induced by pint and pext 


£ — ^int + ^ext • (11) 

The E^ term of £ can be divided into T and L components as 

J J dr[E^^\rf + E^^\rY] . ( 12 ) 

The integral of the L component is evaluated as 

l/drlfid-f = |/drl{£;f'p+p„+f,y 

+ 2Tint ■ Text + 2£^o ^ ■ {Tint + Text}] • (13) 



The terms 


I y dr [ 4 , + 2 ^i„t • £’ext + ■ £i,t] ( 14 ) 

among others influence the motion of internal particles via variational principle. 
Using ([9]), we can evaluate the integral as a sum of the Coulomb potentials 


The integral containing is written as 
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which is zero by using V ■ = 0 and the fact that the upper and lower limits of 

the integral can be taken arbitrarily large outside the isolated matter system. 

Eliminating the variables irrelevant to the motion of internal particles, we obtain 
the Lagrangian for the internal particles 

^ + y dr {-Pint(r)0(r) + Jint(r’) ■ A{r)} 

+Uc + t/ie + I/dr {|^l" - |v X . (18) 

All the particle variables, except for those of Uie, are the ones for internal system. 
Due to the rewriting to make Coulomb potentials explicit, appears only in 
the interaction term. The signs of Uq and are positive here, but there arise 
cancelling contributions from the interaction terms. 

Using = dP^^/dt , we can rewrite the L-fleld related interaction terms as 


/ drl-^),„,(r)^(r) + j!J;’(r)' 

drP,i> . Si"-* + jJ drP,i> . Ad-) , (19) 

where V ■ = —Pint- The first term of the r.h.s. can be further written, by the 

help of eqs.flHl- fTOj) , as 

I drP^^i ■ PT) = _2Uc -U,e +I drP^^i ■ E^^l (20) 

This allows us to rewrite Ci as 

£i = i: + /dr j!J,>(r) • A<’')(r) + / drP<^> ■ E<P - Uc 

The minimum action principle is not affected by the total time derivative term in 
Lagrangian. The possibility of its elimination means that the EM field variables 
do not contain the L field induced by pext- 










Before showing how to eliminate the total time derivatve term, we point out 
the gauge invariance of the many particle Schrodinger equation based on eq. lfTSj) . 
The generalized momenta are 

= TUiVe + eeA{ri) , IIat = , (22) 

and the corresponding Hamiltonian is 

r dA^^^ 

^ ^ + y dr IIat • - 

= Y. l^lPe - ^eA{re)]^ - He - t/ie + [ dr pint(r)0(r) 

£ J 

The last term is the Hamiltonian of transverse EM held (vacuum held) and the 
rest is the matter Hamiltonian in a given EM held 

'^m = Y tt—- eiA{ri)y - He - He + /drp(r)0(r) . (24) 

£ /m£ J 

The Schrodinger equation of the internal particles in a given EM held is 

ih— = (25) 

and its solution allows us to calculate the expectation value of any physical quanti¬ 
ties, especially the induced current density, which plays the role of the constitutive 
equation to determine the EM response together with Maxwell equations. 

It should be noted that this many-body Schrodinger equation is gauge invariant, 
i.e., the gauge transformation {A, 0, T} —)■ {A',0', T'} mediated by an arbitrary 
scalar function x(r, t) 


A' = A + Vx, = T' = exp(i0)T, (26) 

® = ^J dr pint{r)x{r,t) . (27) 


does not change the Schrodinger equation. If we denote the coordinate, velocity, 
and spin of an arbitrary particle as O, we have O^' = exp(i0)(9\[', so that the 
expectation value 

{dj'\fid)m = m{dm ( 28 ) 

of any physical quantity written as a function of (9’s is gauge invariant. 

When we treat materials containing heavy elements, or magnetic species, etc., 
it is required to consider the relativistic corrections, such as spin-orbit interaction 
{'Hso), spin Zeeman interaction (Tfz), etc. For that purpose, we simply add the cor¬ 
responding terms to T-Lq- The spin Zeeman term can be included in the matter-EM 
held interaction term, and the rest in the matter Hamiltonian "Hm- This correction 
does not change the argument about the gauge invariance, since the correction 
terms are written in terms of E and .B [H |2] . 







3 Discussions 


3.1 Special meaning of Coulomb gauge 

In order to eliminate the total time derivative term in the we simply need to 
make a gauge transformation to cancell it, which turns out to be Coulomb gauge, 
as shown below. If we apply the gauge transformation fl26l) to C\, the difference 
5Ci arises only from the interaction terms as 


= / d.|pS . v| + ^ . Vxl = I / d. pS . Vx . (29) 

If we choose y satisfying Vy = —A^^\ the total time derivative term of fl2T]l is 
canceled. The choice Vy = eliminates the L component of A, so that it is 

equivalent to the choice of Coulomb gauge. 

Thus, the Lagrangian in Coulomb gauge is 


£c = E + / drP(L) ■ + J dr ■ A^^^ 


eo 


dr 




and the corresponding Hamiltonian is 


(30) 


Wo = E + /drP<‘-l.pM 
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(31) 


It should be stressed that the EM variables in this Hamiltonian are and 
E^\ which are the gauge independent components of EM held. Combining the 
facts that T-Li gives the gauge invariant Schrodinger equation and that T-Lc is writ¬ 
ten only in terms of gauge independent components of EM held, we directly get 
the gauge invariant EM response from the Schrodinger equation in terms of "He- 
The induced current density is given as a power series expansion with respect 
to these held variables, and their coefficients, i.e., the susceptibilities, are writ¬ 
ten in terms of the eigenvalues and eigenfunctions of the many-body Hamiltonian 
'Hmh = + Uc (including, if necessary, relativistic corrections). The 

merit of this representation is that the susceptibilities are given as separable in¬ 
tegral kernels, which plays an essential role in solving integral equations and also 
carrying out the long wavelength approximation to derive macroscopic constitutive 
equation The arguments given above enforces the foundation of these micro- 

and macroscopic response theories by assuring the gauge invariant nature. 


3.2 Homogeneous L field vs. ’’External” L field 

The L electric feld E^\ the solution of homogeneous Gauss law, has turned out to 
play an essential role in the gauge invariance of many-body Schrodinger equation. 
However, the fact that it has nothing to do with internal and external charge 








densities makes us wonder how it is, or is not, related with the description of 
a matter sample disturbed by an external charge density in cases, such as the 
polarization of a matter sample placed in a condencer, or energy loss spectroscopy 
of an electron beam incident on a matter sample, etc. In these examples, the 
electric held due to charged condenser or electron beam is treated as an external 
electric held inducing a L polarization in the sample. This type of analyses are 
known to work well, but the problem here is how it hts to the gauge invariant 
formalism. 

Obviously, we cannot ascribe the L electric held arising from charged condenser 
or electron beam to We should rather leave as a free L held and regard 

both of the matter sample and the charges inducing the L held as the ’’internal 
system” of the previous sections. In this way we can keep the gauge invariant 
many-body formulation. The gauge to describe this internal system is arbitrary, 
but Coulomb gauge will be of pactical convenience in many cases. Dividing the in¬ 
ternal system into two parts, [A] sample and [B] charged particles on the condenser 
or electron beam, we consider the coupling of the two parts for a given initial con¬ 
dition. In this scheme, we should generally consider the mutual action and reaction 
between [A] and [B]. In the case of condenser, it polarizes the matter sample and 
its polarization changes the capacity of the condenser ahecting the charges accu¬ 
mulated on the condenser. In the case of electron energy loss spectroscopy, electron 
beam induces the matter excitation with L polarization, plasmons for example, and 
the corresponding change occurs in the electron beam, showing the energy and mo¬ 
mentum loss corresponding matter excitation. In the latter case, the energy loss 
function for the beam can be represented by the inverse of L dielectric function 
l/e{k,oj) of the matter 

3.3 Problem of ’’velocity vs. length gauge” 

In the quantum mechanical description of matter - EM held interaction, there have 
been a lot of arguments as to the form of interaction between charged particle and 
EM held, either —E ■ r or {—e/m)p ■ A, since early days They are said to 
correspond to diherent gauges of EM held, and are usually called length gauge and 
velocity gauge, respectively. 

These two forms of interaction are connected by Power-Zienau-Woolley (PZW) 
transformation in the Lagrangian |^, which is specihed by the addition, to the 
Lagrangian in the Coulomb gauge, of the total time derivative of 

F = - J drP{r) ■ A{r) . (32) 

The new Lagrangian equally serves to the quantum mechanical description of the 
system. Through this addition, the interaction term 

Ant = J drJ{r) ■ A{r) , (33) 

becomes 

= J dr{p(T)(r) ■ p(T)(r) + M(r) ■ P(r)} (34) 

via the identity J = {dP/dt) -|- V x M. The interaction term J ■ A stands for the 
’’velocity gauge” and P^'^^ ■ E^^'^ for the ’’length gauge”. 


It is instructive to compare the total time derivative term dF/dt with the similar 
term caused by a gauge transformation, eq. fl2^ . For any choice of y they cannot 
be same, because the T component exists in F but not in the latter. This 
shows that the nonuniqueness of a Lagrangian with respect to the addition of an 
arbitrary total time derivative term is a broader concept than the one due to gauge 
transformation. (This is reasonable from a more general viewpoint, i.e., Lagrangian 
formalism works also in systems which have nothing to do with electromagnetism.) 
For this reason, —E ■ r or {—e/m)p ■ A" is not a problem of gauge. 

4 Summary 

In summary, we have shown that the many-body Schrodinger equation with ex¬ 
plicit Coulomb potential can be given in a gauge invariant form, and that the EM 
response obtained in Coulomb gauge directly gives the gauge invariant one. 
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